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-^ . Abstract. The relative current density n* of "conduction" neutrons in a neutron 

^ ! star crust beyond the neutron drip threshold can be expected to be related to the 

^ I corresponding particle momentum covector pi by a hnear relation of the form n* = 

O ! ^^"^Pj ^^ terms of a physically well defined mobility tensor )C^-^ . This result is 

^ ' describable as an "entrainment" whose effect - wherever the crust lattice is isotropic - 

■^ ' will simply be to change the ordinary neutron mass m to a "macroscopic" effective 

Q_|. mass 772^ such that in terms of the relevant number density n of unconfined neutrons 

we shall have AT*-^ ^ [n/m^j'j^-^ . In a preceding work based on a independent 
particle treatment beyond the Wigner-Seitz approximation, using Bloch type boundary 
j^ I conditions to obtain the distribution of energy Ek and associated group velocity "U^, = 

dSk/dhki as a function of wave vector ki , it was shown that the mobility tensor 

K^ \ would be proportional to a phase space volume integral /C*-^ OC / d A; V^Vf^ 5\Sk ~ 

H I /^} ) where /i is the Fermi energy. Using the approach due to Bogoliubov, it is 

shown here that the effect of BCS pairing with a superfluid energy gap Ap and 



corresponding quasiparticle energy function &k = yi^k ~ /^) + ^f ^^H J^^^ 1^^ ^° 
replace the Dirac distributional integrand by the smoother distribution in the formula 

AT*-^ OC / d k v\v i^, iS.'p / €^ i^ . It is also shown how the pairing condensation gives 
rise to superfluidity in the technical sense of providing (meta) stability against resistive 
perturbations for a current that is not too strong (its momentum Pi must be small 
enough to give 2|j9i'y)j| < €!j^/\£k ~ /^| for all modes). It is concluded that the 
prediction of a very large effective mass enhancement in the middle layers of the star 
crust will not be significantly effected by the pairing mechanism. 



1 Introduction 

A two fluid model of neutron star cores, of the kind that is commonly used in hydrody- 
namical simulations has been developed in the past, assuming a neutron-proton-electron 
composition and using a non relativistic treatment, by the work of many authors PIE]- 
In particular, as explained by Borumand et al. [3], appropriate expressions have been 
obtained for the relevant entrainment coefficients, relating the momentum of the neu- 
tron fluid to the particle current of both neutrons and protons, in terms of the Landau 
parameters in the Fermi liquid theory. More recently, it has been shown by Comer 
et al. j31 E] how the required entrainment coefficients, and the corresponding effective 
masses for the neutrons and the protons, can be evaluated in a relativistic treatment 
- such as is appropriate for deeper layers at densities substantially above the value (of 
about 10^^ g/cm'^) that characterises ordinary nuclear matter - using a relativistic a — ui 
mean field model. The salient conclusion to be drawn from all this work is that in these 
homogeneous fluid layers the effective mass of the neutron will be significantly but not 
enormously reduced (by several tens of percent) below its ordinary bare mass value. 

The present article is part of a newer program of work [7j concerned with the 
previously unstudied problem of evaluating what turns out to be a much more substan- 
tial effective mass modification that is obtained from the corresponding entrainment 
coefficients in the crust at subnuclear densities, above the neutron drip threshold (at 
about 10^^ g/cni'^). In these layers, relativistic corrections are insignificant but the 
issue is complicated by the microscopic inhomogeneity of the medium, in which some 
neutrons can still flow freely but protons are confined to atomic nuclei. As in the outer 
neutron star crust where all nucleons are bound, the nuclei of the inner crust will be 
liquid at high temperatures but will be held by Coulomb forces in a crystalline solid 
lattice in the relatively low temperature range (well below 10^ degrees Kelvin) observed 
in ordinary isolated neutron stars, which according to theoretical considerations (see 
e.g. |H]) should be attained within about a hundred years after the birth of the star. 
Such a low temperature regime will be maintained even in a binary system involving 
accretion provided its rate does not exceed the typical order of magnitude ~ 10^^° 
Mo/yr 0. It has long been generally recognised [TU] that at such low temperatures 
(indeed all the way up to 10^ degrees Kelvin or more ^T]) the neutrons will form a BCS 
type condensate characterised by a superfluid energy gap, and - for the reasons dis- 
cussed in the penultimate section of this article - will therefore be able to flow through 
the lattice for a macroscopically long time jT2j , without resistive or viscous dissipation. 
The evaluation of the entrainment coefficients for such a superfluid flow is of particu- 
lar astrophysical interest, because relative motion of the effectively free "conduction" 
neutrons through the inner crust lattice is believed to be an essential element in the 
mechanism responsible for observed pulsar frequency "glitches" . 

In order to obtain the quantitative information that is needed, a ffist step has been 
the development by the present authors of a microscopic derivation [U] of an appropriate 
two fluid model for the inner crust regime, on the basis of a simplified non-relativistic 
description of the underlying nuclear physics in which (as in the previously cited work 



j2l El on the higher density regimes pertaining to deeper layers of the star) thermal 
corrections and the effects of superfluid pairing were neglected. Recent numerical work 
[7j has shown that this treatment implies an enormous enhancement by the entrainment 
of the effective mass of the free neutrons in the middle layers of the inner crust (a result 
that contrasts with the rather moderate diminution predicted for the effective neutron 
mass in the deeper core layers). The purpose of the present work is to evaluate the 
adjustments - which turn out to be small - that will result from allowance for the 
superfluid pairing of the neutrons. 

Until now, quantum theoretical analysis of neutron superfluidity has mainly con- 
centrated on static configurations (in either infinite medium ^3] or inhomogeneous 
systems |14| ITK| ITT)]), meaning states for which no current is actually flowing relative 
to crust. Even at densities substantially beyond the neutron drip threshold it should 
still be possible to obtain a reasonably accurate description for the static case by using 
the so called Wigner-Seitz approximation that treats the neighbourhood of each ionic 
nucleus as if it were isolated in a sphere whose diameter is determined by the nearest 
neighbour distance. However for the treatment of more general - stationary but non 
static - configurations involving relative conduction currents it is absolutely necessary 
to use a more realistic description in which the artificial Wigner-Seitz type boundary 
conditions are replaced by the natural Bloch type periodicity conditions that would be 
desirable for higher accuracy even in the static case. 

The use of appropriate Bloch type periodicity conditions is routine in terrestrial 
solid state physics [T7|, but has so far been applied to neutron star matter only in a 
simplified zero temperature independent particle treatment [HIIZ] (of the kind applicable 
well below the Fermi temperature, but most appropriate for a young neutrons star that 
has not yet fallen below the critical temperature for the onset of superfluidity) for which 
the neutrons are considered to move as independent particles without allowance for 
the pairing interactions responsible for the superfluid energy gap that (in cool mature 
neutron stars) allows the currents to persist. 

A simplified treatment of this kind has been used to show that the middle layers of 

a neutron star crust will be characterised by a very low value for the relevant mobility 

tensor in the formula n^ = JC^-^Pj for the current Ti^ ^ TiV^ of unbound neutrons 

(which will be present above the "drip" density of the order of 10^^ g/cm'^) with number 

density n , mean velocity V^ and momentum per neutron Pi . Throughout this paper 

summation is understood over repeated covariant and contravariant coordinate latin 

indices, for instance JC^-^Pj = J2j JC^-^Pj- In the independent particle treatment, the 

mobility tensor was shown |^6^ to be given by a volume integral over the space of Bloch 

momentum covectors ki that is expressible in terms of a Dirac distribution with support 

confined to the Fermi surface - where the relevant energy function S^a with a band 



index a, is equal to the chemical potential ^U - in the form 

J^'' = 7^ E / vlAa S{Ska - /i} d'k , (1) 



{27ry 



a 



(in which, as throughout this work, we use braces - as distinct from ordinary brackets 
- for functional dependence, in order to avoid possibly confusion with simple multipli- 
cation) where the relevant group velocity distribution is given by the usual formula 

The main purpose of this article is to show how the preceding independent particle 
treatment can be generalised to allow for BCS type pairing using an approach of the 
kind pioneered by Bogoliubov. Since the relevant temperature range (substantially 
below 10^ degrees Kelvin) is very small compared with the critical value (of the order of 
10^ degrees Kelvin) for the pairing condensation it will be justifiable for us to continue 
to use the zero temperature limit, T = , in which thermal correction effects are 
entirely ignored. One of the main motivations for this work is to check the robustness 
of the conclusions obtained from the simple treatment described above, particularly the 
prediction of a very low value for the mobility tensor, which is interpretable as meaning 
that the corresponding effective mass m^ = n/3/C\ will become very large compared 
with the ordinary neutron mass. 

Our conclusion is that as a first step towards a more accurate treatment, in cases 

for which the superfluid pairing can be characterised just by a gap parameter Ap 

the relevant integral over the Fermi surface will need to be replaced by a phase space 

volume integral given in terms of the quasi-particle energy 



€ka = V(^fca - f^y + ^F (3) 

by the new formula 

JC'^ = 7A-,T.JvL<^d'k, (4) 

in which the expression for the group velocity v\^ is the same as in the absence of 
the pairing gap. It is the diminution of this group velocity that is responsible for the 
enhancement of the effective mass, which on average should therefore not be greatly 
affected by the phase space smearing effect produced by the superfluid pairing. This 
comes from the fact that the expectation value in the superfluid phase of one particle 
quantities, namely the particle current density in the present work, are not very different 



from their normal (non superfluid) value as shown by Leggett ^H] in the context of 
superfluid helium 3. 

Although the effect of the pairing is not so important for the evaluation of the ef- 
fective mass, it is of course important for the property of superfluidity itself. Much 
of the contemporary literature on the underlying mechanisms for "superconductivity" 
in astrophysical contexts deals only with purely static configurations, in which the es- 
sential question concerns the existence of a condensed state characterised by a finite 
energy gap. A secondary purpose of this article is to go back to the question of super- 
conductivity in the strict technical sense, which refers not to static configurations but 
to stationary configurations involving the relevant motion of a current of some kind - 
electric in ordinary laboratory metals, but neutronic in the case of interest here. The 
essential issue is that of the (meta) stability of such a current against (small) pertur- 
bations of the kind that in the "normal" case would produce resistive damping. The 
original defining property of a superconductor is that it should be able to support a cur- 
rent that will be characterised by such metastability provided it does not exceed some 
finite critical value, beyond which "normal" dissipation will of course set in. It will be 
confirmed in the penultimate section of this article that within the framework of the 
simple theoretical model used for the present work this condition of superconductivity 
will indeed be satisfied, with a critical maximum current value that is estimated to be 
safely large compared with what is required for the relevant applications to astrophysical 
phenomena such as pulsar glitches. 

2 Hamiltonian for the independent particle limit 

The idea is to start on the basis of a second quantised formalism in terms of local 
fermionic field annihilation and creation operators ip and ip^ depending on space 
position coordinates X^ in a unit volume sample, and on a spin variable G taking 
values j and | subject to anticommutation rules of the usual form 

[^a{r}, V^a'{r'}]+ = , [^V{r}, ^V{r'}]+ = , (5) 



[^Ur},^.'{r'}]+ = ^..'^{r,r'}, (6) 

using a quadratic Hamiltonian operator of the form 

/\ /\ ,^ /■y. /\. y^ 

ind iiit ' ind kin pot ' \ / 

in which the interaction term H .^^^ will be absent in the independent particle limit 
corresponding to the kind of model used ^ in our preceding first quantised treatment. 



In this independent particle limit, only the kinetic and potential contributions are 
present and will be assumed to be given (neglecting possible spin dependence for sim- 
plicity) by integrals over the unit volume sample under consideration of the form 

^,,„ = / dV H, Jr} , 7Y, Jr} = E ^^.{r}H,>4r} , (8) 

^pot = /dVT^^Jr} , TY^Jr} = l^{r}E^V{r}^4r} , (9) 

where l~Ly^^^ is a self adjoint differential operator in the category specified in terms of 
a gauge covector flj by an expression of the familiar form 

Ha = -Y'iV, + la,) \ , (V, + la,) (10) 

2m® {r} 

in which 'j'^-^ is the Euclidean space metric and 7?2®{r| is interpretable as a micro- 
copic effective mass, which is usually found to have smaller values inside crustal nuclei 
|14j . The covector with components tti is a gauge field allowing for the possibility of 
adjustment of the phases of the field operators ^/^cril"} • Iii applications to particles 
with non zero electric charge ( e say) such as the electrons in an ordinary terrestrial 
superconductor or the protons in the deeper layers of a neutron star, the presence of 
such a field (taking the form flj = eAi ) would be necessary for the treatment of 
magnetic effects, but in the uncharged case of the crust neutrons with which we are 
concerned here it will always be possible to work in the standard gauge for which this 
covector is simply set to zero, Cij ^ , which means that we simply take 

K.=no. (11) 



"kin 



The potential V^{r} and the microscopic effective mass 777-® |r| (as those deduced 
from contact two body interactions of the Skyrme type in the Hartree Fock approxi- 
mation) are supposed to represent the averaged effect on the neutrons of the nuclear 
medium and in particular of the ionic lattice. A periodic crystalline type lattice will 
be assumed, which implies that these functions should be invariant with respect to any 
lattice translation vectors : 

y{r + rej = F{r}, m®{r + TeJ = m®{r} (12) 

for any triad of integers i°' (a ^ 1, 2, 3 ) in which the lattice basis vectors e^ 
may be interpreted as representing the interionic spacing in the solid case that will be 

6 



relevant at very low temperature, but should in principle be taken to be much larger (so 
as to generate a giant cell interpretable as a typical mesoscopic average over a locally 
disordered configuration) for applications above the relevant melting temperature, at 
which it is to be expected that (unlike the weaker electron pairing mechanism in or- 
dinary terrestrial superconductors) the superfiuid neutron pairing mechanism will still 
be intact. 

The mass function ?T2®{r} in the specification of the kinetic contribution to the 

independent Hamiltonian will also be involved in the specification of the corresponding 

neutron current density operators, which will be given, for each value of the spin variable 

a , by 

K{^] = ^-||^7^^' (^l{r}V,V^,{r} - (V,V^l)V^,{r}) . (13) 

One of the main objectives of the present work is to obtain a practical way of evaluating 
the mean value of the total current, as given by the space averaged operator 

^' = E^:, #:=/dVn;, (14) 



as a function of the associated momentum in a stationary state that is non static (and 
therefore non isotropic, since the mean current will characterise a preferred direction) 
but uniform of a mesoscopic volume, meaning one that is large compared with the 
interionic spacing but small compared with the macroscopic lengthscales characterising 
the star crust thickness or even the intervortex separation. 

It is to be noted for future reference that this current can be used to express the 

adjustment that will be required in cases when it turns out to be more convenient to 

work with the gauge adjusted operator l~La rather than 1~Lq in the kinetic contribution 

flllj) : it can be seen that this kinetic contribution will be given in the small a limit by 

^M„ + «^^* = E /dV V^t^{r}7i„V^,{r} + 0{\a\^] , (15) 

subject to the usual assumption that we are using periodic boundary conditions to get 
rid of a boundary term produced by an integration by parts using Green's theorem. 

3 Representation by Bloch states 

Subject to the usual Bloch type boundary conditions for a mesoscopic material sample 
of parallelopiped form - with a unit volume that is taken to be very large compared with 
the elementary lattice cells under consideration - the independent particle Hamiltonian 



will determine a complete orthonormal set of single particle states (^fcQ,{r| , labelled 
by a wave covector kj taking discrete values on a fine mesh inside the first Brillouin 
zone and a band index a , satisfying the Floquet-Bloch theorem ^Tj 

ifka{r} = u,a{r}e"'-\ (16) 

using the abbreviation A; • r = /c^ X* , where Ukai^j satisfies the ordinary lattice 
periodicity conditions 

Uka{r + rea} = Uka{r} . (17) 

These wave functions are normalised as follows (using to indicate complex conjuga- 
tion): 

J d'riflArJn^ir} = 6MSap . (18) 

Setting the Bloch wave vector ki in place of tti in the definition (jlOj) the eigenvalue 
equation can be usefully rewritten in terms of the ordinarily periodic functions Uka as 

{Hk + V)Uka{r} = SkaUka{r} . (19) 

From the spin independence of the potential (by neglecting spin-orbit coupling terms), 
the phases can be chosen in such a way that we shall have 

"Pkai^} = ^-ka{r} , ul^{r} = W_fca{r} . (20) 

These Bloch states may be employed in the usual way as a basis for the specification 
of corresponding position independent annihilation and creation operators, C^afc and 
C aka ) subject to anticommutation relations of the standard form 

[Caka, Ca'll3]+ = , [c^aka, C^a'lf3]+ = , (21) 



Caka, Ca'lf3]+ = ^aa'hl^afS , (22) 



in terms of which the original position dependent annihilation and creation operators 
will be given by 

V^a{r} = E ^ka{r}Caka , '^^ajr} = E ^ka{^}c'aka ■ (23) 

k,a k,a 



It is evident just from the ortho normality conditions (J18p that the spin dependent 
number density operator defined by 

na{r} = ^V{r}^4r} . (24) 

will have a space integral 

#, = /dVn,{r}, (25) 

(over the unit volume sample under consideration) that will be given by 

k,a 
It can similarly be seen from the defining conditions ()23p that integrated value of the 
independent particle contribution, 

H,^, = /dV?Y^„,{r} , 7^^„,{r} = H,Jr} + H^Jr} , (27) 

which is interpretable in the absence of the interaction contribution as the total energy 
operator, will be expressible in standard form as 

-^ind = H ^ak flaka , (28) 

where ^o-fea is the Bloch wave vector dependent particle number density operator 
given by (j^ . 

To get an analogous formula for the mean current (over the unit volume sample 

under consideration) as given by the operator (J14j) , we take its expectation value 

(|iT;,|>=/d='r(K{r}|), (29) 

for a state | ) satisfying the simplicity condition that except for the diagonal con- 
tributions characterised by O" = (7, /^ = ki and a ^ (3 the contributions of 
the expectation values \\c^akaCa'li3\) will vanish - or be negligible to the order of 
approximation under consideration - it can be seen that we shall obtain the formula 

m)=T.{\naka\)vl^. (30) 

k,a 

in which the relevant velocity will be given by 

4. = / dV ,^r^lj-,l'K^la{^Wm.a{v] - ^.a{r}V,^L{r}) . (31) 

This expression ()31|) for the velocity vector v\.^ can easily be shown to be mathe- 
matically equivalent to the well known, albeit less intuitively obvious, group velocity 
formula that is given in terms of the single particle energy introduced in ()19|) by (j21) • 



4 Characterisation of conducting reference state 

The (zero temperature) states in which we are interested are those that minimise the 
expected total energy (|-^|/ subject not only to the usual constraint that there should 
be a fixed given value of the corresponding total expected particle number 

(l^l) = E(l^.|), (32) 

a 

but also, since we are concerned with non-static - conducting - stationary configura- 
tions, to the requirement that there should also be a fixed given value of the expected 
total 

(i*i> = E(|s;d , (33) 

a 

of the current defined by (fT^ . 

Imposing these constraints by the introduction of corresponding Lagrange multipli- 
ers /i and Pi , the problem will effectively be that of unconstrained minimisation of 
the combination 

(|ff'(„)|) = (|ff|>-„(|#|)-p,(|fti>, (34) 

in which we introduce the notation 

H[py = H' - p,fi\ H' = H-im. (35) 

In the absence of the pair coupling term H .^^^ , the quantity to be minimised reduces 
to the form (|//|^^|^^||) with 

^LM = ^L-l^^^^ K. = H..-m. (36) 

It can be seen from (|TT| and (|T5|l that, for a small value. 

Pi = Hi , (37) 

of the momentum, its effect will be given to first order in the magnitude \q\ of the 
corresponding wave number covector just by substituting the gauge adjusted operator 
l~L-q in place of 1~Lq in the relevant differential formulae. Thus, in particular, it can be 
seen that the appropriate modification of the eigenvalue equation ([11^1 for the required 
replacements S UAka ^ind UipXi^^y of 8ka and Uka "will be just the substitution of 
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Tik-q for T~('k ; which evidently means that to this first order of accuracy we shall 
have 

^{p}ka = £{k-q)a , (38) 

and 

U{p}ka = U(^k-q)a ■ (39) 

One thereby obtains the formula 

^{pjkai^} = e"^"'(p^k-q)a{Y} , (40) 

for the corresponding modification of the single particle states (|T6|l , which in turn, 
by the analogue of ()23p , determine correspondingly adjusted annihilation and creation 
operators, Crp|^y^,Q, and c' {pi^^.Q, , in terms of which (PUJ) can be rewritten in the 
equivalent form 

'^cr 2^ '^{p}aka 5 '^{p}aka C [pjaka^{p}aka • \'^^) 

k,a 

It can thus be seen that, as the analogue of ()28|) , the effective Hamiltonian H r \, 
will be given by the formula 

^Mind = 5Z ^'{p}ka n{p}ka , (42) 

a,k,a 

with 

^{p}ka = ^{k-q)a ' (43) 

to the linear order accuracy in p with which we are working. At this order, it can be 
seen from (j21) that (jlSj) may be rewritten as 

^'{p}ka = ^'ka - Pi '^ka , ^'ka = ^ ka " f^ ■ (44) 



The expectation value of the quantity given by (j42|l will evidently be minimised by 
a reference state vector | ) ^ I {/(,«}) that is chosen (as a function of the multipliers fl 
and J9j ) in such a way that the expectation \{j_i^p}\n^p'^g-ka\{^,p}) has its maximum 
value, namely 1, whenever ^|„i;tQ, is negative, and its minimum value, namely zero, 
whenever ^ rpU,Q, is positive. It can thus be seen from (jlHj) that the effect of the current 
will consist just of a uniform shift of the distribution in momentum space by an amount 
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given by the infinitesimal momentum covector Pi . Such a state is characterised by the 
conditions 

n{p}aka\{^i,p}) = n{p}aka\{^i,p}) ■ (45) 

with the eigenvalues given as a Heaviside distribution by 

n{p}aka = ^{-^'{p}ka}- (46) 

It can be seen that this state |{u,p}) satisfies the condition for applicability of the 
analogue of (ISUj) , and hence that the expected mean current value 

K = {{,,pM\h.,p}) (47) 

will be given (in accordance with our previous evaluation |B] in a first quantised frame- 
work) by 

K = T.vLH-^'{p}ka}- (48) 

fc,a 
In the linearised weak current limit with which we are working, it can be seen from 
(j21) and ()44|) that this will be expressible to first order in terms of the static limit value, 

riaka = t^i/U - ^ka} , (49) 

of the distribution PUJ) as 



5 Bogoliubov treatment of pairing 

Up to this point what has been done in the present article is just to translate the work 
of our preceding article [Hj from first quantised to second quantised formalism. The 
motivation for this translation is that a second quantised treatment is indispensible 
for the next step, which is to go beyond the independent particle model used in the 
preceding work by including allowance for pairing interactions. 

In the inner crust in which we are concerned here, dripped neutrons are expected to 

be paired in spin singlet states ^Uj, as the usual Cooper pairs of electrons in terrestrial 

superconductors. A standard way of allowing for pairing interaction in a mean field 

model (assuming for simplicity a contact two body interaction as it is the case for 
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conventional superconductor ^\ and a common practice in nuclear physics [20]) is 
thus to take the interaction contribution in (0) to have the form 

TY^Jr} = A{r}^\{r}V^\{r} + A*{r}^|{r}^^{r} , (51) 

where A|r} is a position dependent complex potential that, in a "self consis- 
tent" model, should be expressible in terms of the abnormal density expectation value 

(l^/^llrj^/^ilr}!) in the relevant reference state |) . 

The mean complex phase of the function A{r} is subject to an indeterminacy 
that can be resolved by fixing the phase in the specification of the wave operators. In a 
static configuration one would expect that this coupling potential A{r} would share 
the ordinary lattice periodicity property fll2|l and moreover that the phase should be 
adjustable in such a way as to ensure that A becomes real. 

Instead of using the representation ()lfj|l in terms of the simple Bloch wave functions 
y^fcail"} ) ill the approach introduced by Bogoliubov one seeks a more general repre- 
sentation whereby the single component Bloch waves are replaced by two component 
Bloch functions with components V'fcal^} ^^^ V^fcal^J ^^^^ are characterised with 
respect to corresponding ordinarily periodic functions W^.^{r} , and M^^j,{r} by 

These functions are used for replacing the original representation ()16p by a mixed 
particle-hole representation involving a new set of position independent quasi-particle 
annihilation and creation operators ^^rk and 7 V/fc in terms of which we shall have 

i^rU] = E (^L{r}7Tfca - ¥^fca{r}7iL) . (53) 

k,a 

and 

fc ^ ^ 

where the new operators satisfy anticommutation relations of the standard form 

[Tafca, 7a'lf3]+ = , [tW, l^<7'll3]+ = , (55) 



[l^aka, ^a'll3]+ = ^aa'hl^aP ■ (56) 

As a result, consistency with (0) and © entails the relations: 

[V^Hr}, M^'}U = n^. r'} = E ^ti^WhU'} + <{r'WL{r} , (57) 

k 
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[i;\{r}, V^t^{r'}]+ = = E ^LUMaW} " AUMaW} ■ (58) 

k 

The purpose of such a Bogohubov ansatz is to enable us to choose the new functions 
V-'fcal^} ^^'^ V^fcal^J ^^ ^^"-^^ ^ ^^y ^^ ^° simphfy the expression for the total effective 
Hamiltonian, which will be given for a static configuration by 

H' = Jd^rH'ir}, (59) 

with 

n'{r} = E 4{^}KM^} + A{r}^|{r}^|{r} + A*{r}V^^{r}^^{r} , (60) 

(7 

in which the independent particle contribution is given by 

K. = ^in. - /^ , (61) 

where, as before, /i is a Lagrange multiplier, whose purpose when we apply the vari- 
ation principle, is to impose the constraint that the expectation of the total integrated 
number density should be held fixed. It is to be remarked that in the presence of the 
pairing interaction term, the number operator n will no longer exactly commute with 
the Hamiltonian, which implies that the state that minimises the expectation of the 
effective Hamiltonian obtained in this way will not be an exact eigenstate either of the 
particle number or of the energy. 

The simplification of the Hamiltonian ()6U|) can be achieved by taking the functions 

(/9^Q,{r} and V?fcQ,{r} to be solutions of the coupled set of differential equations (known 

as the Bogoliubov-de Gennes equations in the condensed matter field |2I]) given by 

■"'■4) (:;:)--(£)■ 

in which the eigenvalue €^ka is what will be seen to be interpretable as the rele- 
vant quasiparticle energy. This system can be written more explicitly in terms of the 
ordinarily periodic functions ul.^\r} , and 'W^Q,{r} introduced in (jK^ as 

using the notation of (J6T|l , where 

V' = V-fi. (64) 
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The foregoing specification is incomplete, because the condition of satisfying ()(J2|1 will 
evidently be preserved by interchanges of the form 

^ka ^ ^-ka . €!ka ^ "^'/ta , (65) 

but this ambiguity is resolved by adoption of the usual postulate that the eigenvalues 
be positive, 

6^fca > . (66) 

To fix the normalisation of the solutions, which will automatically satisfy the integral 
relations expressible - restoring the explicit reference to the position dependence - as 

/ dV ^L{r}^?/,{r} = / dV ^L{r}^/V{r} , (67) 

the amplitude of the (automatically mutually orthogonal) solutions is chosen so that 
/ dV«{r}^?^{r} + ^tirWifAr}) = hi^^p • (68) 

The foregoing ansatz eliminates all the undesirable terms, reducing the effective 
Hamiltonian operator to the remarkably simple form 

H' = J2 €^ka {l^akalaka " siu^Oka) = J2 ^ ka {^OS^Oka " 1 akalKka) , 
a,k,a a,ka 

(69) 
in which Oka is the relevant Bogoliubov angle, as defined, for each value of the 
wavenumber covector k^ and band index a by 

cos^et„ = / dV <{r}^L{r} , sin^et„ = / dV vi;{r}^L{r} ■ 

(70) 



By minimisation of the expectation of the operator ()69|) one obtains the required 
condensate reference state | ) = |{//,|/ , which is characterised by the condition 

7afca|{M}) =0, (71) 

expressing absence of all the quasiparticles created by the operators 7 Vfco; • 

The quasiparticle operators can be written in terms of the particle operators re- 
membering equation ()23p and the orthonormality condition (fTHj) as 

C]ka = Y. [Uka.lplW ~ ^ka,ll3l Up) , (72) 

1,(3 ^ ^ 
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where we have introduced the matrices 

t/w/3 = /dVv^*,Jr}vp?^{r}, (74) 

V^W/3 = /dV^%Jr}^g{r}, (75) 

which from the properties of Bloch wave functions reduce to 

Uka,l[3 = 5kl Uka,kl3 , Vka,W = ^-kl Vka-kp ■ (76) 

It is to be noted that ifka i '^%a ^^"^ V^A-a ^^^ ^^ Bloch wave functions associated 
with the same Bloch wave vector (hence having the same phase shift whenever trans- 
lated from one cell to another) but are solutions of different equations. Inserting these 
expressions into the number density operator TT-o-fca introduced in ()26|) , it is readily 
verified that its expectation value in the superfluid ground state is given by 

{{pL}\naka\{pL}) =^\Vka,-kl3\ ■ (77) 

p 

Remembering that '^ka are the single particle states of the independent Hamilto- 
nian (jF)T|) , with energies £^^.q, , it can be seen from the Bogoliubov equations (jH^ that 
the expression (fTTjl is equivalent to 

13 yp-kp + o^.^) 

where 

AL,^/3 = / dV ^L{r} A{r}^g{r} . (79) 



6 The BCS ansatz 

Since (particularly for the middle layers of a neutron star crust, where the effective 
mass enhancement is likely jHl El to be most important) we are still far from having a 
sufficient knowledge of the solutions (/9^.Q,{r} for the independent particle model, it 
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will evidently take some time before we can hope to obtain a complete evaluation of 
the solutions for the coupled equations for '^^^1^} ^^"^ '/'fcal^} using an accurate 
estimate of the coupling coefficient A{r} . In the meanwhile, as an immediately 
available approximation, offering the best that can be hoped for as a provisional estimate 
in the short run, we can use an ansatz of the standard BCS kind, which means adopting 
the prescription 

Uka,lf3 = COS 9ka5kl5a/3 , Vka.ip = Slli 9 ka^-k^l^ap ■ (80) 

Comparing with (J53|l and (J54|) , the Bogoliubov particle-hole doublet reduces to 

^laiA = fka{^} COS Oka , ^kaiA = ^kai^} sin Oka , (81) 

where the single component wave functions Lpka\^j are the (more easily obtainable) 
independent particle eigenf unctions, which can be seen from the preceding work to be 
specifiable as solutions of the simple Schroedinger type equation 

Knd^ka = S'ka^ka , (82) 

where, in the static case under consideration at this stage, we simply have 

^'ka = ^ka - ^ (83) 

where Eka is the ordinary Bloch energy value as introduced in (fTII|) . 

It can be seen that the ansatz ()81|) will provide an exact solution in the limit for 
which the relevant coupling field matrix elements 

^ka.ip = J dV ^L{r} A{r}^,,{r} , (84) 

reduce to diagonal form, so that we have 

^ka,lf3 = ^ka^kl^aP , (85) 

using the notation 

Afca = ^ka,ka ■ (86) 

The relation ()8(i|l will be a good approximation when A/^.q, remains close to a fixed 
value Ap (which can be taken without loss of generality to be real and positive 
by choosing the relevant phase) in the neighbourhood of the Fermi surface, and it 
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will evidently hold exactly when the coupling constant is uniform, so that A{r} = 
Ap = Afca-. In the general case, as a result of the periodicity of A{r} the pairing 
field matrix elements will automatically be diagonal in phase space, namely A^q, ;^ = 
Ski/S.kQ.^kj3. However the pairing interactions may couple single particle states belonging 
to different bands and it will only be an approximation to neglect those contributions 
when, for instance, A|r} is a field of the radially dependent form that has been 
obtained jI3] within the Wigner Seitz approximation. Actually the only non vanishing 
matrix elements are those relating independent single particles states belonging to the 
same irreducible representation of the space group 'TZ , which means that only the band 
states having the same symmetry properties may be coupled. Subject to the validity of 
fl85p , the BCS ansatz (jHT|l will reduce the Bogoliubov system of differential equations 
(Ifi2j) to a purely algebraic eigenvalue system whose solutions have the well-known form 



Ska = 'JS'tJ + AL , (87) 

cos^(),„ = %^ , sin^fl,-„ = %^ . (88) 

It can be seen that the ansatz ()8H) has the effect of reducing the Bogoliubov trans- 
formation to the simple form given by 

C-aka = COS 9ka 1 -aka + C^ sin Oka iJ-ka ' (89) 

which is equivalent to taking 

J -aka = COS Oka C-aka " ^^ sin Oka cj-ka ■ (90) 

It follows from this that for the state \{n}) characterised by ()7H) , the expectation 

values of the Bloch wave vector dependent number density operators flaka introduced 

in (pF)|) will be given by 

{{^l}\naka\{^l}) = siia^Oka , (91) 

This result is interpretable as expressing the effect commonly described as a smearing 
of the Fermi surface, whereby the smoothed out Bloch wave vector space distribution 
(1911) replaces the hard cut off expressed by the Heaviside formula ()49|1 that applies in 
limit when the pairing interaction is ignored. 



7 Formula for the mobility tensor 

When the static contribution characterised by (j61|) is extended by the inclusion of the 
current constraint term proportional to the momentum covector Pi = Tiqi in the 
effective energy (jH^ , it can be seen from ()15|) that as in the independent particle limit, 
its effect to first order will be entirely taken into account by merely making the gauge 
adjustment ai = —qi, in the kinetic energy operator ()1()|1 , which means changing ki 
to ki — Qi in equation (J63|l . The first order effect of the current will therefore be 
given, according to (jiij) , by the adjustment 



^ka '^ ^{p}ka — ^{k-q)a 5 (92) 

for the single particle energy, and by the ensuing set of infinitesimal transformations 

7aka ^ 7{p}aka , |{m}) ^ \{fi,p}) ' (^3) 

of quasi-particle operators and state vector, while particularly, in the framework of the 
BCS approximation based on the neglect of interband couplings, the Bogoliubov angles 
introduced in (jHTj) will undergo a corresponding adjustment 

Oka i-^ 0{pyka ■ (94) 

As in the absence of pairing, in the strict BCS case characterised by a fixed gap value, 
A/jQ, = Ap, the result will still be describable just as a uniform displacement 5ki = Qi 
in the space of Bloch wavevectors ki . 

As the adjusted version of ()H()|1 , it can be seen that for any state | ) satisfying the 

simplicity condition that except for the diagonal contributions characterised by (7 ^ 

(7, li = ki and Q; = (3 the contributions of the expectation values \\& {p\akaC-{p\a'lj3\) 

will vanish - or be negligible to the order of approximation under consideration - the 

mean current defined by (jTTjl will be given for each spin value by the formula 

K = T.vla{\^{p}crka\). (95) 

k,a 

In the framework of the BCS approximation this formula will be applicable, in partic- 
ular, to the conducting reference state | ) = |{/ijj}/ , so that by the adjusted analogue 
of fj89|) the ensuing replacement of the formula ()48|) , for the mean current in this state, 
will be obtainable from the substitution 

{{Pi,p}\'^{p}cjka\{p,p}) = sin (9{p}fca , (96) 
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which leads to the expression 

n' = T.K = ^T.vL^^%p}ka. (97) 

f^ k,a 

for the corresponding total current. 

Since the total current evidently cancels out in the unperturbed static state \{^i\) , 
the quantity given by (jHTj) will be expressible to first order, in the weak current limit 
with which we are working, as 

n-Z}^Vf.^Pj . (y«j 

k,a Opj 

The conclusion to be drawn from this is that the value of the current will be given to 
linear order by an expression of the same general form 

n' = PjKL'^ , (99) 

as in the absence of pairing, but with required mobility tensor now given by an expres- 
sion of the form 

^■- = 2E4 /'"f"""^ (100) 

It follows from ()92|) that in this small \p\ limit we shall have 

-4ao (101) 



^^{p}ka _ ^^{p}ka _ dSka 



dpi Thdki hdki 

and hence that the partial derivative in (j^JHj) can be evaluated in the BCS approximation 



as 



in which sin Oj^^a is given as a function of the quantity £j^^ = Ska ~ f^ cind of Ap 
by ()88p . The mobility tensor will therefore be expressible as 

^'^ = -2E ^y'"^ «. (103) 



(104) 



in which, by ()88p , the relevant coefficient will be given by 

oSka '^€^ka 
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The translation of the discrete summation formula ()103|) into the language of con- 
tinuous integration (in the limit in which the size of the mesoscopic cell is much larger 
than the lattice spacing) is given by (jlj) . 

Except near the base of the neutron star crust where the nuclei may acquire exotic 

(e.g. "spaghetti" or "lasagna" type) configurations, it is to be expected that the mobility 

tensor will have the isotropic form 

IC'^ = JCY^ , (105) 

where 

p- 1 i^ij 2 (9 (sin ^/^-n) 2 2 i 1 /in^\ 

/C = 37^J■/C ^ = -- L —Q^ Vka , ^ka = 7*J^fca^L , (106) 

It is to be observed that subject to the BCS approximation of uniform coupling, 
meaning that there is a constant gap parameter, ^ka ^ ^F ? the formulae p03|) and 
f)lU6|) will be convertible, using integration by parts, to the form 

rii _ O \p ^^^ ^^-Q' ^ ^ka ^ _ 2 ^ sin Oka d Ska /-.^^N 

This latter formula is useful for the evaluation of the corresponding effective mass 
m^ as defined by 

m* = n//C, (108) 

in terms of the relevant total particle number density as given by the prescription 

^ = E(l^a|)=2EsinX.«, (109) 

<^ k,a 

in which, if we only wish to count unbound neutrons, the summation should be taken 
only for values above a lower cut off below which the states are bound so that the 
corresponding values of the velocity Vka will vanish. 

The concept of an effective mass has traditionally been a source of confusion as 

different definitions have been used in different contexts. Moreover in solid state physics 

one is often more interested in electric charge (not mass) whose transport is related to 

the electric field Ej via an Ohm type law as 

f = en' = a'%, (110) 
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where Ti* is the electron current density, and (7*-^ is the relevant electric conductivity 
tensor. While this conductivity tensor cr*-' has the advantage of relating macroscopic 
measurable quantities, it depends on the dynamical evolution of the medium unlike the 
newly introduced mobility tensor lO'-^ , on which the effective mass 772^ is defined. The 
electric conductivity tensor will be given by an expression of the form 

a'J = e^rJC'^ (111) 

in which T is a timescale characterising the rate of decay (by various scattering pro- 
cesses) towards the zero current state that is the only locally stable configuration in the 
"normal" case. The kind of superconducting case with which we are concerned may be 
described as a limit in which the relevant timescale r is infinite, so that the conductiv- 
ity (7*-' will also be infinite, even though the mobility tensor JC^-^ has a well behaved 
finite value as in the "normal" case. However it is important to understand that the 
reason why the relevant timescale T is effectively infinite in the superconducting case 
is not because relevant scattering cross sections are small (as in the case of a "normal" 
good conductor) but rather because the current carrying configuration is locally stable 
(with respect to scattering processes that may be quite strong) in a superconducting 
state, for reasons that will be reviewed in the next section. 

It is to be remarked that the formula for the mobility tensor p()H|l is very similar to 

the formula obtained without pairing correlations, the Heaviside unit step distribution 

being merely smeared. In particular the same velocities appear in these formulae. 

One might have naively guessed that apart from the particle state distribution which 

is smoothed, the relevant velocity would have been given not by the ordinary group 

velocity "y^^^, given as the momentum space gradient of the energy distribution Ska 

by (0) but by the analogously defined quantity vl.^ obtained by substituting €!ka iii 

place of Ska 5 namely 

Actually this latter "pseudovelocity" is interpretable as a mean velocity between parti- 
cles and holes, since €^ka is the energy of a quasiparticle which is a mixture of particles 
and holes. More specifically, when (as in the simple BCS case for an homogeneous sys- 
tem) the gap parameter is independent of the momentum, this modified velocity will 
be given by the expression 

^L = ^L^, (113) 

from which it can be seen that v]_.^ will vanish at the Fermi surface characterised by 
^ka = /^ ; where the number of particles is equal to the number of holes. 
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In the limit for which, in so far as the unbound neutrons are concerned, the effect 

of the crustal nuclei is small (either because the nuclei occupy only a small part of 

the volume, as will be the case just above the neutron drip transition, or because the 

nuclear surface is very diffuse, as will be the case near the base of the crust) we shall 

have 

1 d'^^ka ^ 1 ,^^^. 

U^ dkidkj m® ^'^ ' 

where Ul® is the uniform mass scale appearing in the kinetic energy operator, which 
will be comparable with, but for precision somewhat less that, the ordinary neutron 
mass m . It can be seen by comparing (jlU7p and (jlU9|) that in this approximately 
uniform limit the effective mass for the unbound neutrons will be given simply by 
m^ = TTl® , regardless of whatever the value of the gap parameter A may be. 

8 Superconductivity property and critical current 

An unsatisfactory feature of the rather profuse contemporary literature dealing with 
various kinds of what is commonly referred to as "superconductivity" in astrophysi- 
cally relevant contexts (including such exotic varieties as colour superconductivity in 
quark condensates) is the rarity of any serious theoretical consideration of the actual 
property of superconductivity in the technical sense, meaning the possibility of having a 
relatively moving current that is effectively stable, or in stricter terminology metastable, 
with respect to small perturbations - such as would normally give rise to a dissipative 
damping mechanism of a resistive or viscous kind. 

In the astrophysical literature concerned with pulsars it has generally been taken 
for granted that neutron currents of the kind considered in the present work actually 
are characterised by superconductivity in the sense of being metastable with respect to 
relevant kinds of perturbation. In this section we shall investigate the conditions under 
which this supposition of metastability is indeed valid. The issue is that of the stability, 
for small but finite values of the momentum covector Pi , of the superconducting refer- 
ence state I ) = \{np}) that is characterised by the minimisation of the combination 



The conducting state |{u,p}) was derived by minimising the energy expectation 
(|i7|) subject to the condition that the particle number expectation \\n\) and the 
current expectation (|7I*|) were held fixed. It is physically reasonable to suppose 
the particle number expectation (|^|) really will be preserved under the conditions 
of chemical equilibrium that are envisaged in the relevant applications, but no such 
consideration applies to the current expectation (|7i*|) which in a "normal" state 
would tend to be damped down by many conceivable kinds of scattering process. The 
physically pertinent question is therefore whether Ijuj?}) will still minimise \|-^|) 
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with respect to small relevant perturbations - subject of course to the preservation 
of the particle number expectation (|^|/ as before - when the prior assumption of 
preservation of (|?T*|) is abandonned. Subject to the particle number conservation 
condition 

(5(1^1) = 0, (115) 

this stability requirement is equivalent to the condition of minimisation of {\H |) 
meaning that any admissible perturbation must satisfy 

6{\H'\)>0, (116) 

where according to the notation introduced in (jH^j) , 

H' = H[^^+P,n\ (117) 

According to the reasonning of the previous section, the relevant adjustment of 
(J7U|) will give us 

-^M = 51 €! {p}ka[T{p}akal{p}aka- sin 0{p}ka) , (US) 

SO that the specifications ()117|) and (j^^j) provide the variation formula 

^(1^1)= E {€^{p}kJ{\j\p}akal{p}aka\)+PtVlJ{\n{p}aka\)) , (119) 

cr,fc,Q- 

in which, for the BCS case, it can be seen from ()87|1 that we shall have 



€{p}ka = JSl,,^^ + Al. (120) 



'{p}ka - \j^{p}k 

In this BCS case, the action on the conducting state | {;/,»}) of a typical quasiparticle 
creation operator 7 (pjn.Q- will provide only three non vanishing terms in the sum 
(I119|l , namely those given by 

5{\l\p}]kal{p}^ka\) = 1, (121) 

together with the number variation contributions 

^(I^MT-tal) = COS^%}fca , 5{\fl^p}^_ka\) = -Sm^e{p}_ka ■ (122) 

It follows from the symmetry properties v\,^ = —v\^^ and Oka = ^-fca that 
the explicit dependence on the Bogoliubov angle will cancel out at first order in the 
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net energy contribution provided in (J119j) by such an excitation, so that this energy 
contribution will be positive if and only if 

€{p}ka+P^vl^>0. (123) 

It is to be noted that such an individual quasiparticle excitation will in general violate 
the requirement (jlisp to the effect that the number of real particles should be conserved, 
but it is evident from ()122|) that such violations may have either sign and so can be 
cancelled out by the combined effect of two or more elementary excitations. What, in 
a stable case, can not be cancelled out is the combined effect of several quasiparticle 
energy contributions in ()119|) : it can be seen that the quasiparticle energy contributions 
will always add up to give the positive result needed for stability provided the inequality 
(I12H|I is satisfied for all admissible modes. 

The stability condition ()123|1 that we have derived for a BCS type inhomogeneous 

superconductor is consistent with Landau's classical treatment based on vaguer heuristic 

arguments in the context of superfluid Helium 4 [21j. Since the quasiparticle energy 

€' {plka is always positive, it is clear that the stability condition (jl23|) will always be 

satisfied ii Pi vl.^ > 0. Therefore an instability in the superfluid conducting state can 

only occur when J9j "U^^-^ < 0. Rewriting the inequality ()123p as €^ [pika > ~Pi '^\a^ 

squaring both sides and substituting the expression (HHj) for ^|p|fcQ, in ()12U|) we see 

that in the BCS case there is a remarkable simplification (which does not seem to have 

been pointed out before) whereby the terms that are non linearly dependent on the 

momentum covector Pi cancel out, so that the superfluid conducting state will be 

stable if the following inequality is satisfied 

2j.,i;L4a<^L- (124) 

This can be rewritten in terms of the "pseudovelocity" introduced in ()112|) as 

p,v\^<\€ka. (125) 

which simplifies to the following form whenever the BCS gap is independent of the 
momentum 

P^^X\ri{€L])<n. (126) 

The inequality ()124j) is evidently verified for elementary excitations above the Fermi 
level for which S^,^ > 0. In the derivation of the inequality ()124|1 we have assumed 
that PiV\^ < 0. Actually there will always be some modes for which this is satisfied 
since whenever PiV\^ > we shall have p^'L'LfcQ, < . 
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The stability condition (jl23|) (which has to be satisfied for all modes) can therefore 
be restated as the requirement that the magnitude p of the mean particle momentum 
covector Pi lies below some critical threshold J9c 

P<Pc (127) 

in which, for an approximately isotropic distribution depending only on the magnitude 
k of the wavenumber covector k^ , the critical value Pc will be given by 



Pc-r k:^ l4J + a^ . (128) 




Since £^^,q, vanishes on the Fermi surface, it is clear from ()128|) that Pc will also 
vanish - so that there will be no phenomenon of superconductivity - not only when the 
gap A^Q, vanishes everywhere, but even when it vanishes just in the neighbourhood 
of the Fermi surface. When the mean value Ap of the gap at the Fermi surface is 
non-zero but small compared with the other relevant energy scales - as will typically be 
the case - it can be seen that the minimum in ()128|1 will be attained for energy values 
differing from the Fermi value by a small but finite positive or negative amount that will 
be given approximately by ^^,q, ~ zt Ap . In such a case, it follows that the critical 
momentum value ()128|) will be expressible in terms of the mean value "Up of the group 
velocity magnitude V^a ^^ the Fermi surface by the approximation 

Ap 
Pc « — • 129 

Introducing the critical velocity as Vq = P(,/ni^ , the criterion p29j) can be written 
as 

Ap 
v^ « — — . (130) 

In the limit of an homogeneous medium for which m^ = 771® , the critical velocity 
reduces to an expression which is commonly found in the literature concerning homo- 
geneous electron superconductivity in metals [SHj, namely 

^(0) 

^(«) « -f , (131) 

(using the superscript ("^ to indicate what would be obtained for uniform values of the 
microcoscopic effective mass 777-® and potential V ) where kp is the radius of the Fermi 
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sphere. It must be emphasized however that the critical velocity of an inhomogeneous 
superfluid (such as the neutron superfluid in the inner layers of a neutron star crust) 
will differ from the estimate ()13H) by a factor 



Vr Af St. 



v'^' A';' 4"^ ' 



(132) 



where op and op are the Fermi surface areas in the inhomogeneous and homogeneous 
cases respectively. Since the opening of band gaps in the single particle energy spectrum 
decreases the Fermi surface area, the critical velocity is therefore expected to be smaller 

than the expression ()13H) assuming that Ap ~ Ap . 

For a gap of the order of an MeV, in a region where the kinetic contribution to 
the Fermi energy has a typical value of the order of a few tens of MeV, the formula 
fjl29|) implies a critical momentum value corresponding to a kinetic energy of relative 
motion of the order of hundreds of keV per neutron. This is comparable with the total 
kinetic energy of rotation in the most rapidly rotating pulsars. However the relative 
rotation speeds of the neutron currents that are believed to be involved in pulsar glitch 
phenomena are very much smaller - by factors of 10"'^ or even far less [23 ^ than 
the absolute rotation speeds of the neutron star. In all such cases it may therefore 
be concluded that the superfluidity criterion ()129j) will be satisfied within an enormous 
confidence margin. 

It is to be remarked that in the more thoroughly investigated context of laboratory 
superfiuidity |21] Landau's simple linear formulation of the stability problem in terms 
just of phonons provides only an upper limit on the critical momentum whose true 
value is considerably reduced by the less mathematically tractable - since essentially 
non linear - effect of rotons. Analogous considerations presumably apply in the present 
context. This means that although our present treatment places the estimate ()129p on 
a sounder footing than was provided by any previous work of which we are aware, it 
should still be considered just as an upper bound on the true critical value which is 
likely to be substantially reduced by non linear effects whose mathematical treatment 
is beyond the scope of the methods used here. Despite this caveat, the prediction 
of genuine superconductivity in the context of glitching neutron star crusts should be 
considered to be very robust. The justification for such confidence is that - according 
to the considerations outlined in the preceding paragraph - the relevant magnitudes 
of the neutron currents in question correspond to values of the neutron momentum p 
that are extremely small compared with the order of magnitude given by ()129|1 . For 
such very low amplitude currents there is no obvious reason to doubt the validity of 
conclusions - including estimates of effective masses, as well as the prediction of genuine 
superfluidity - that are based on the simple kind of linearised treatment used here. 
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9 Conclusions 

In the middle layers of the crust, where the effect of inhomogeneities will be important, 
our previous analysis neglecting the effect of the superfluid gap lead to the prediction [H] 
that there will be a strong "entrainment" effect whereby the value of TTl^ will become 
very large compared with m . This prediction has now been confirmed by an analysis 
[7j based on the phenomenological model of Oyamatsu et al. [211 where values as large 
as m^jm ^^ 15 has been found at a baryon density 77-5 ^ 0.03 fm . Our present 
analysis indicates that this conclusion will not be significantly affected by taking the 
relevant pairing gap A into account. Thus although the pairing is essential for the 
actual phenomenon of superconductivity, on the other hand, in so far as the effective 
mass is concerned, neglect of the effect of superconductivity pairing will indeed be 
justifiable as a robust first approximation, at least for moderate values of the gap 
parameter compared to the kinetic contribution to the Fermi energy. 

The unimportance of pairing from the point of view of entrainment, which has been 
usually assumed (see for instance Borumand et al. j3]) and explicitely shown in the 
present work, can be seen from the consideration that, when Aq-^ is small compared to 
the Fermi energy /i , the coefficient ()l()4j] will be very small everywhere except in a thin 
layer with width of the order of ^ak near the Fermi surface locus where S^k ^ f^i 
which means that when the coupling is weak its effect will be entirely negligible. In 
sensitive cases for which the geometry of the energy contours near the Fermi surface is 
complicated by band effects, a moderately strong pairing effect might make a significant 
difference by smoothing out variations of the mobility tensor as a function of density, 
but does not seem likely that this smearing effect would make much difference to the 
large scale average properties of the mobility tensor. In other words the effective mass 
is expected to be much more sensitive to band gaps than to the pairing gap. The main 
reason is that A appears only in the number density distribution whereas band gaps 
(resulting from Bragg scattering of dripped neutrons by crustal nuclei) have a strong 
influence upon the neutron velocity "y^^ which is vanishingly small in this case. 
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